Impurity Scattering in Luttinger Liquid with Electron-Phonon Coupling 
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We study the influence of electron-phonon coupling on electron transport through a Luttinger 
liquid with an embedded weak scatterer or weak link. We derive the renormalization group (RG) 
equations which indicate that the directions of RG flows can change upon varying either the relative 
strength of the electron-electron and electron-phonon coupling or the ratio of Fermi to sound veloc- 
ities. This results in the rich phase diagram with up to three fixed points: an unstable one with a 
finite value of conductance and two stable ones, corresponding to an ideal metal or insulator. 
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Interacting electrons in one dimension are known to 
form a Luttinger liquid (LL) characterized by power- 
law correlation functions [1 . This characteristic feature 
of the LL has been established via conductance mea- 
surements and a scanning tunneling microscopy (STM) 
both in carbon nanotubes [2] and semiconductor quan- 
tum wires |3j. Embedding a potential impurity into the 
LL leads [H [5] to a universal (i.e. impurity-independent) 
power-law (in temperature T) suppression of the trans- 
mission amplitude through the LL and suppression of 
the tunneling density of states (TDoS) near the impu- 
rity, with the latter fading away with the distance [5J [7] . 

The electron-phonon (el-ph) coupling in addition to 
the (Coulomb) electron-electron (el-el) repulsion in the 
LL is known to result in the formation of two polaron 
branches with different propagation velocities (see, for 
example, [8]). In the present Letter we show that embed- 
ding a single scatterer into such an el-ph liquid results in 
a rather rich phase diagram: depending on the relative 
strength of the el-el and el-ph coupling and on the ratio of 
the Fermi to sound velocities, the system can be an ideal 
metal (with conductance g = 1 in the units of e 2 /h as 
in the absence of impurities [9] ) or an ideal insulator for 
any impurity strength, or to be in an intermediate state 
from which it can flow cither to the metallic or insulating 
limit, depending on the impurity strength. 

To show this we consider both the weak back-scatterer 
(WS) and weak link (WL) limit, the latter correspond- 
ing to a weak tunneling coupling between the two halves 
of the LL. In the phononless case the WS and WL am- 
plitudes are known to scale at low energies as e 7 ~ and 
e 1+ , respectively For the repulsive el-el interaction the 
exponents 7_ = K — 1 < and 7+ = K~ l — 1 > 
where K < 1 is the Luttinger parameter. Therefore, 
the renormalization group (RG) flows always go in the 
direction of stronger scattering (weaker transmission) re- 
sulting at T — in the insulating phase for any scatterer 
[3]. The el-ph coupling leads to j± — > j±, with each of 
the exponents 7± changing sign at different values of bulk 
parameters. In the WS limit which has been previously 
considered [ID] , the change of sign of 7_ was reversing 
the RG flows indicating that the weak scatterer becomes 



irrelevant for the sufficiently strong el-ph interaction and 
the LL thus remains in the metallic state. But the con- 
siderations of both the WS and WL limits presented here 
show that the change of signs in the exponents j± hap- 
pens at different values of the parameters, resulting in 
the rich phase diagram described in detail later. 

To get the results, we employ the functional bosoniza- 
tion formalism in form developed in [7J. This allows us 
first to include the el-ph interaction which leads to elec- 
trons dressing with phonons, i.e. the formation of po- 
larons, and only after that to bosonize the action. Prior 
to considering the embedded scatterer we describe the 
polaron formation in a pure LL, reproducing the known 
results [B] in form convenient for further considerations. 

We consider a model of ID acoustic phonons linearly 
coupled to the electron density. The phonon spectrum is 
assumed to be linear with a cutoff at the Debye frequency, 
u D = cq D (a straightforward modification for the 3D or 
optical phonons will be described elsewhere) . Integrating 
out the phonon field in the standard way results in sub- 
stituting the dynamical coupling, V(£) = Vo + T>(£), for 
the screened Coulomb interaction Vq in the LL action: 

Sll=E UMO^MO-l [dZV(0n 2 (0. (1) 
v =±i J AJ 

Here d v = d t + T]v p d x , the (spinless) electron field is de- 
coupled into the sum of left- (rj = — 1) and right- (77 = 1) 
moving terms, -0(0 = V'R^) 6 '^ + tph(S,) e ~ lPFX with 
£ = (x, t), and n = (V'rV'r + V'lV'l)- We use the Keldysh 
formalism [11 implying the time integration along the 
Keldysh contour in Eq. |l]) and below. The free phonon 
propagator 2?(£) in V(£) = Vo + 25(0: can be defined by 
the Fourier transform of its retarded component, 
— 1 2 

V o( w > ?) = — 2 2~ ' W 9= C <?' UJ + = LJ + lO, (2) 

where c is the sound velocity, a p h is the dimensionless el- 
ph coupling constant, and vq — (ttVf) is the free spin- 
less electron DoS. The form of the LL action in Eq. 
implies neglecting the electronic backscattering. This re- 
mains justified in the presence of the el-ph coupling fol- 
low temperatures, T <C w D . 
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The next step is the Hubbard-Stratonovich transfor- 
mation which decouples the n 2 term in the action and 
results in the mixed fermionic-bosonic action in terms of 
the auxiliary bosonic field ip minimally coupled to ip: 

SeS = ~l j^vV^ip + ij^id^-ip)^. (3) 

We gauge out the coupling term by the transformation 



corresponding to el-el repulsion (which becomes stronger 
with the el-ph coupling) and attraction. It is the ex- 
istence of these two modes that leads to a rich phase 
diagram when a scatterer is embedded. 

Following Kane and Fisher [1] we consider two types 
of scatterers (assuming them pinned at x — and not in- 
volved in lattice vibrations): a weak backscatterer (WS) 
or a weak tunneling link (WL). The WS action is 



V(0=¥>(0. (4) S ws = A /d^ R (t)VL(t)+c.c, m = ^ = 0,t), (7) 



The Jacobian of this transformation results [7] in substi- 
tuting V' 1 + II for V- 1 in Eq. where II = II R + n L 
is the one-loop electronic polarization operator (exact for 
the LL [12]). with 11^ (£) = ig v (£)g v (—Q and the free elec- 
tron Green function defined via the Fourier transform of 
its retarded component as gZ(u), q) = [u) + — rjv F q]^ 1 . 
The phase 9 V is related to the auxiliary field 99 by [13] 



where is the bosonic Green functions that resolves 
. Its retarded component coincides with the free 
fermionic (while the Keldysh components are natu- 
rally different). 

The Green function of the interacting polarons is not 
gauge-invariant with respect to the transformation Q 
and depends on the correlation function lU^i = (9 v 6 ri r): 

The retarded Fourier component of £/„,„' is found as 



UL'(u,q) = 



J++r)'v F q V {uj 2 + - uj 2 ) + U 1 



a ph UJ 2 



, (5) 



m > u + -r)v F q (uj\ - v 2 + q 2 ) (lj 2 + - v 2 _q 2 ) 

where v± are velocities of the composite bosonic modes: 



v 2 + c 2 ± \J (v 2 — c 2 ) 2 + Aa v 2 c 2 



(6) 



Here v is the speed of plasmonic excitations in the 
phononless LL, v = v F (l + vqVo) 1 ^ 2 = v F K~ x , where K 
is the standard Luttinger parameter. Without phonons 
(c = 0), one has w_ = and v + = v, so that in this case 
(as well as for uj > uj d or a p h = 0), U m ' reduces to the 
usual LL plasmonic propagator. 

We assume the parameter a = a p ^K 2 in Eq. Q obey- 
ing the inequality a < 1 to avoid the Wentzel-Bardeen 
instability [14 corresponding to v 2 _ < (with the thresh- 
old shifted from a p ^ — 1 for a pure el-ph model to 
a ph = l/K 2 > 1). Then the velocities v T of the slow 
and fast composite bosonic modes in Eq. (|6| obey the 
inequalities u_ < c, v < v + . These modes mean that the 
LL in the presence of the el-ph coupling becomes two- 
component, with the effective Luttinger parameters 



fast 



v F /v + < K < 1 



K f 



low 



v F /v- > 1, 



with Ao being a bare backscattering amplitude. 

The WL action contains the tunneling term linking the 
two halves of the LL, labeled by 1 and 2: 



t dt ipxfa 



Ipa = V'aR + Ipal 



(8) 



with to being a bare tunneling amplitude. We begin with 
the WS case, noticing in advance that the WL case is dual 
to it as well as for the standard LL (4]. 

The gauge transform Q replaces Ao in Eq. ^ with 

A(t) = A e i ^W- e «W] ) 9 v (t) = 9 v (x=0,t). (9) 

Integrating out all x ^ fields results in the quadratic in 
6(t) action with the scattering term (|9| as in the phonon- 
less problem [3]. The difference is that the correlation 
function is governed by the composite bosonic modes, 
Eq. (pj), and its retarded Fourier component is 



VM>' 



dq 
2tt 



i / ^U m i{u;q) 



, (10) 



where the correlation matrix 7^ is found from the 
straightforward integration above as follows: 



7 LL = 7r R = a(K fast ) £+ + "(^siow) f- , 
7 LR = 7 L r = a(-Kfest) £+ + a{K slow ) . 



(11) 



Here £± = \ [l ± (v 2 - c 2 ) (v 2 + - 
functions of the two modes and 



are the weight 



a(K) = UK- 1 + K) - 1, a{K) = \{K^ - K). (12) 

In the absence of phonons £+ — 1 and £_ — and the fast 
mode contribution reproduces the results of [4] . 

The only change in the x = action which is due to 
the el-ph coupling is the appearance of the dimensionless 



matrix 7^^/ in Eq. ( 10 ). Thus the RG analysis for the WS 
case is similar to that in [4] but the RG equation for the 
backscattering amplitude A acquires a different prefactor 



defined by Eq. ( 11 ), 7_ = 7 LL - 7 



di X(E) 



(1 - K)X(E) , E>lo d 
-7-A(£), E<uj d 



(13a) 



Here E is a running cutoff and I = \uEq/E with Eq ~ e F 
being the bandwidth. 



(a) 



(b) 



FIG. 1. Weak tunneling (7+) and weak scattering (7-) exponents as functions of the Luttinger parameter K < 1 for a weak 
(a) or strong (b) el-ph coupling, (a p h ^ 1, respectively). The dotted lines represent y± for the phononless LL, the solid and 
dashed lines correspond to /3 F = v F /c = 0.5 and /3 F — 1.5, respectively. The latter is shown only for the strong coupling where 

for /3 F > [a p h(a p h — 1)] ^ 2 the exponent 7+ does not change sign within the stability region, K < a^ 1 ' 2 . 



The WL case can be treated with the same 



: ac- 
tion defined by Eq. (10 1 (rather than going to the dual 
</>-action [3] which is inconvenient for our mixed repre- 
sentation PS])- The gauge transform Q changes the 
tunneling action (^sj) by replacing t/j a —> i/v^e" 79 " and 
t Q -> to&^- t) : where 0± = \{Q aK - 9 aL ). The RG 
analysis of the tunneling term yields the RG equation 



dit(D) 



(l-K-i)t(E) 
|13a|), with 7+ = 7 LL 



E>u T . 
E <w E 



(13b) 



dual to Eq. 

The RG equations ( 13 ) should be solved with the initial 
conditions X(E = E ) = A and t(E = E ) = t Q . The 
solutions for e < w D are given by 



He) 
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En 



m 

to 



Eo 



7+ 



7+ 



(14) 



The exponents j± can be rewritten (with /? = v/c) as 
7ct = KfjK a - 1 , a = ±1 , (15a) 



1 



(15b) 



In 



while 7± are obtained with k± = 1 in Eq. ( |15a 
the a bse nce of the el-ph coupling we have "f a — j a and 
Eqs. (14) reduce to the standard ones, A(e) ~ (e/E'o)^" 1 
and t(e) ~ (e/E Q ) K _1 , i.e. the backscattering ampli- 
tude increases (the WS case) and the tunneling ampli- 
tude decreases (the WL case) with e — > for any K < 1. 
This shows that embedding an arbitrary scatterer results 
in the LL becoming (at T = 0) an ideal insulator. 

The el-ph coupling, however weak, changes drastically 
the above conclusion leading to the possibility of a metal- 
insulator transition with changing K. 



Indeed, it follows from Eqs. ( 13 1— ( 15 ) that RG flows 



can change directions depending on the values of K, a p h 
and /3 F = v F /c. It is easy to verify that k + (K) < 1 while 
K - {K) > 1 which results in j± changing sign (see Fig. 1) 
when the interaction strength K equals K ± , with 

7+ (^+)=0, 7_(X_) = 0; K_<K+<1. (16) 

Then for K + < K < 1 both j± change sign so that 
the WL amplitude t(e) increases while the WS ampli- 
tude A(e) decreases with e — > in Eqs. (14), indicating 
[T3] that adding a single scatterer does not change the 
metallic nature of the pure LL (where the dimensionless 
conductance g = 1) - in contrast to the phononless case. 

Decreasing K (i.e. increasing the el-el interaction), we 
enter the region < K < K + , when the WL ampli- 
tude still decreases - as in the phononless case, but also 
the WS amplitude decreases - opposite to the phonon- 
less case. This indicates that a strong scatterer (i.e. weak 
link) results in the insulating behavior (the 3 = fixed 
point) while a weak scatterer leaves the LL in the metal- 
lic phase (the g = 1 fixed point). Thus there should 
exist an intermediate unstable fixed point with a finite g, 
separating the metallic and insulating regimes. 

Finally, for a strong enough el-el interaction, K < 
the RG flows in Eqs. (14) remain for both the WL and 



WS amplitudes qualitatively the same as for the phonon- 
less case, so that any scatterer makes the LL insulating. 
Explicitly, we find K± for a p h <C 1 as follows: 



l-K n 



1 



c 



This means that in this limit there exist all the three 
regimes described above, see Fig. 1(a). 

When the el-ph coupling is not weak (which is the 
case for carbon nanotubes [16 ) our considerations are 
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A =15 




FIG. 2. Phase diagram: the dotted line is the boundary of 
the stability area, a p h = K~ 2 , the solid and dashed lines cor- 
respond to K^:, respectively. The light-grey area represents 
the insulator (both y± have the same sign as in the absence 
of the el-ph coupling). In the dark-grey area (where only 
7_ changes sign) there exists a line of unstable fixed points 
separating the metallic state for a weak scatterer from the in- 
sulating phase for a strong one. The hatched area (reachable 
only for a p h < a* h , as illustrated by some ttp h ) corresponds 
to the purely metallic phase where both y± change signs. 



— 1/2 

bound by the stability requirement, K < a h so that 



for ctph > 1 the Luttinger parameter is confined to the 
region K < a ph J 
pendence on /3 F 



1/2 



In this case there is an essential de- 



It is easy to see from Eqs. (15) and 



(161 that K+ > a I 12 for > [a ph (a ph - 1)] so 



that for such /3 F the purely metallic regime (with both 
7± changing sign) is no longer accessible, as illustrated 
by the dashed lines in Fig. lb. On the phase diagram 
(Fig. 2) all the three regimes described above exist only 
for a P h <a* ph =\ + yi + 4/ ft. 

There are other properties of the LL strongly affected 
by the el-ph coupling. First, 7 + is also the edge expo- 
nent of the TDoS at the boundary (edge) of the wire, 
I'cdgc ~ £ 7+ , the TDoS also experiences, depending on 
the parameters, a transition from vanishing to divergent 
with e — > 0. Then, transport properties of the LL with 
resonant or antiresonant impurity [17], or of the disor- 
dered LL [IS], or of the LL out of equilibrium Q35] also 
experience qualitative changes with allowing for the el-ph 
interaction, as we will show elsewhere. 

To summarize, we have shown that the el-ph coupling 
qualitatively changes the phase diagram of the LL with 
a single impurity. The change is not reducible to a re- 
definition of the Luttinger parameter K: the existence 
of slow and fast polaron modes with different weights in 
different regimes results in different RG flows for a weak 
scatter and a weak (tunneling) link. The resulting phase 
diagram (Fig. 2) has, depending on the parameters of 
the problem, regimes corresponding to purely metallic or 
purely insulating behavior and an intermediate regime 
with two stable fixed points (ideal metal and insulator) 



and one unstable, finite-conductance fixed point. 

This work was supported by the EPSRC Grant 
T23725/01. I.V.Y. and I.V.L. gratefully acknowledge 
kind hospitality at the Abdus Salam ICTP. 



[1] S. Tomonaga, Prog. Theor. Phys., 5, 544 (1950); J. M. 
Luttinger, J. Math. Phys., 4, 1154 (1963); F. D. M. 
Haldane, J. Phys. C, 14, 2585 (1981); J. von Delft and 
H. Schoeller, Ann. Phys., 7, 225 (1998). 

[2] Bockrath et at, Nature, 397, 598 (1999 ); Z. Yao et al, 
ibid., 4 02, 273 (19 99); H. Ishii et al, \ibid., 426, 540| 
(2003)| J. Lee et al., |Phys. Rev. Lett., 93, 166403 (2004)f 

[3j O. A uslaender et al., Science, 295, 825 (20 02); E. Levy 
et fflZ., |Phys. Rev. Lett., 97, 196802 (2006)| 

[4] C. L. Kane and M. P. A. Fisher, Phys~Rev . Lett., 68, 
1220 (1992); |Phys. Rev. B, 46, 15233 (1992)| 

[5] K. A. Matveev, D. Yue, and L. I. Glazman, Phys. Rev. 
Lett., 71, 3351 (1993). 

[6] S. Eggert, |Phys. Rev. Lett., 84, 4413 (2000)] 

[7] A. Grishin, I. V. Yurkevich, and I. V. Lerner, Phys. Rev. 
B, 69, 165108 (2004); I. V. Lerner and I. V. Yurkevich, 
in Nanophysics: coherence and transport, Les Houches 
Summer School Series, Vol. 81, edited by H. Bouchiat 
et al. (Elsevier, New Yo rk, 2005) p. 109. 

[8] D. Loss and T. Martin, [Phys. Rev. B, 50, 12160 (1994)| 
L. Mathey, D.- W. Wang, W. Hofstetter, M. D. Lu kin, 
and E. Demler, |Phys. Rev. Lett., 93, 120404 (2004)] 

[9] D. L. Maslov and M. Stone, Phys. Rev. B, 52, R5539 
(1995). 

[10] P. San-Jose, F. Gu inea, and T. Martin, |Phys. Rev."!^] 
I 72, 165427 (2005)| 

[11] J. Rammer and H. Smith, Rev. Mod. Phys., 58, 323 



(1986); A. Kamenev and A. Levchenko, Adv. Phys., 58, 
197 (2009) 

[12J I. E. Dzyaloshinskii and A. I. Larkin, Zh. Eksp. Teor. 
Phys., 65, 411 (1973). 

[13] After the gauge transformation Q the remaining 
fermionic part can be bosonized with ip ~ e lrf> . This would 
allow one to rewrite the action in the standard, fully 
bosonized form in terms of = 6 + <j> and cj>. However, 
we will use the mixed fermion-boson representation. 

[14] G. Wentzel, | Phys. Rev., 83, 168 (195"Tj| J. Bardeen, Rev. 
Mod. Phys., 23, 261 (1951). 

[15] When either t(e) or A(e) increases as e — > 0, the one-loop 
RG gives only an indication of behavior in this limit. 

[16] A. De Martino and R. Egger, |Phys. Rev. B, 67, 235418| 
(2003)| S. Sapmaz, P. Jarillo-Herrero, Y. M. Blanter, 



and H. S. J. van der Zant, New J. Phys., 7, 243 (2005) 
[17] C. L. Kane and M. P. A. Fisher, Phys. Rev. B, 46, 
7268(R) (1992); D. G. Polyakov and I. V. Gornyi, 68, 
035421 (2003); A. Furusaki and K. A. Matveev, Phys. 



Rev. Lett., 88, 226404 (2002); Y. V. Nazarov and L. I. 
Glazman, [9l7l2 6804 (200 3)] I. V. L erner, V. I. Yudson, 
and I. V. Yurkevich, 100, 256805 (2008) , M. Goldstein 

and R. Berkovits, |104, 106403 (2010)1 

[18] T. Gia marchi and H. J. Schulz, |Phys. Rev. B, 37, 325| 
(1988)| 



[19] D. B. Gutman, Y. Gefen, and A. D. Mirlin, [Phys. Rev.| 



Lett., 101, 126802 (2008); Phys. Rev. B, 80, 045106 
(2009); 81, 085436 (2010). 



